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Abstract 

I introduce several simplified schemes for the approximation of the self-consistency condition of 
the dynamical cluster approximation. The applicability of the schemes is tested numerically using 
the fluctuation-exchange approximation as a cluster solver for the Hubbard model. Thermody- 
namic properties are found to be practically indistinguishable from those computed using the full 
self-consistent scheme in all cases where the non-interacting partial density of states is replaced by 
simplified analytic forms with matching 1st and 2nd moments. Green functions are also compared 
and found to be in close agreement, and the density of states computed using Pade approximant 
analytic continuation shows that dynamical properties can also be approximated effectively. Ex- 
tensions to two-particle properties and multiple bands are discussed. Simplified approaches to the 
dynamical cluster approximation should lead to new analytic solutions of the Hubbard and other 
models. 

PACS numbers: 71.10.-w, 71.27.+a 
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I. INTRODUCTION 



A major effort in condensed- matter theory goes towards the development of new tech- 
niques for the study of correlated-electron systems, such as the dynamical mean-field theory 
(DMFT). When studying the Hubbard model, DMFT is often implemented with a simple 
analytic form of the self-consistency condition that corresponds to using a Gaussian, semi- 
circular or Lorentzian non- interacting density of states This has led to analytic work 
investigating, for example the Mott transition |31]. Results from dynamical mean-field theory 
are expected to be applicable to 3D systems, but since DMFT only considers local processes, 
physics in low dimensions cannot accurately be investigated. 

There have been several approaches that extend the core ideas of DMFT so that low- 
dimensional systems can be studied, including the dynamical-cluster approximation (DC A) 
which systematically reintroduces spatial fluctuations into the dynamical mean-field theory 
^. DCA is an extremely powerful method which has been applied to a wide range of model 
systems jsj], but it can be computationally expensive, and intensive numerical simulation 
may lack intuition. Carrying out analytic work with the dynamical-cluster approximation 
is difficult owing to the coarse-graining step which completes the self-consistency, which 
consists of a partial sum over momenta in small regions of the Brillouin zone. The partial 
sum can not be carried out analytically. 

By neglecting momentum conservation at the interaction vertices, DMFT replaces the 
self-energy, T,{k,u), by its momentum-independent counterpart. In 2D, the spatial fluctu- 
ations which were neglected in the DMFT are expected to contribute significantly. DCA 
includes these fluctuations by dividing up the Brillouin zone into sub-zones which obey the 

n n 

lattice symmetry |J] or more recently using Bett's clusters |6|]. This coarse graining step 
corresponds to a mapping onto a cluster-impurity model. Within each of these sub-zones, 
the self-energy is assumed to be momentum independent, so the coarse graining of the Green 
function can be written 

G(Ki,z) = Y i . (1) 

Here A; is a vector that lies within the sub- zone as shown in Fig. [TJ The sum over momenta 
is normally carried out numerically. 

Since the self-energy is momentum independent within each sub-zone, it is possible to 
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FIG. 1: Schematic of the dynamical cluster approximation. The Brillouin zone is split into sub- 
zones, within which the self-energy is assumed to be constant. The centers of the zones are at the 
momenta Ki and a vector k is defined which relates the center of a sub-zone to a location in that 
sub-zone. 



replace the sum over momenta with a partial density of states, leading to the relation ^, 

DKX^)de 



(2) 



where DkX^) is the non-interacting fermion partial density of states (DOS) for the sub-zone 
centered about (see figure [2]( a)). 

In both DMFT and DCA, the self-energy and the coarse-grained Green function are also 
related through a modified Dyson equation, 



g,\K,, z) - G-\K,, z) = T.{K.,, z) . 



(3) 



Tom QQ{Ki,z), which may be 
Cluster impurity problems 



The self-consistency is closed by calculating the self-energy 
interpreted as the host Green function of an impurity model 
have been extensively studied, therefore a number of different approximations for the self- 
energy are readily available. 

Here, I suggest a simple analytical form for the dynamical-cluster self-consistency step, 
Eqn. [21 I use results from the Hubbard model Tl] to establish if the simplified forms lead to 
acceptable results. The Hubbard model is one of the simplest non-trivial models of electronic 
correlations in the solid state, and examines a single band of electrons, hopping between 



3 



lattice sites with amplitude t, and interacting via a Coulomb repulsion, U, represented by a 
Hamiltonian, 



(q^ creates an electron on site i with spin a and rijo- is the number operator). This model can 
only be solved exactly in one dimension Q| , but has been extensively examined in the infinite- 
dimensional limit corresponding to DMFT j^, and in a wide range of other approximations, 
including DCA [5|. 

This paper is set out as follows: In section |lll I suggest and summarize several approx- 
imate forms for the self-consistency step of the DCA equations. I test the accuracy of the 
new schemes in section IIIII Possible extension for two-particle properties is discussed in 
section [IV] and for multiple bands in section |Vl A brief discussion can be found in section 



II. ANALYTICAL APPROXIMATIONS TO DCA 

An appealing aspect of DMFT is its potential for analytic work. Three forms for the non- 
interacting density of states are commonly used: (1) a Gaussian, which corresponds to a 
hyper-cubic lattice with high dimensionality, (2) a semicircular density of states which relates 
to a Bethe lattice with large coordination number and (3) a Lorentzian, which decouples 
the self-consistent equations linking DMFT with the Anderson impurity model. In the case 
of the DCA, the non-interacting partial density of states (partial DOS) has a complicated 
form, which prevents easy anal5d;ic work. In the very large cluster limit each partial DOS 
tends towards a 5-function form. Thus in principle, any approximate partial DOS with the 
properties of a 5 function has the correct large cluster behavior (for example one might 
replace the partial DOS with a Gaussian). One may go further and match the mean (e^J 
and variance (tj = a/ e^Ki — (ej<rj^) for the approximate and exact partial DOS, to get the 
correct large cluster asymptotic behavior (since the low order moments dominate the form 
of the self-consistent condition in that case). Here the bar indicates an average over the 
momentum states in the subzone centered at Ki. Smaller cluster sizes are likely to be most 
useful for analytic work of the type carried out with DMFT, thus it is of interest to determine 
whether such a procedure is useful when Nc = 4, which is the smallest cluster size which 
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TABLE I: Offsets and widths of tlie partial DOS for a 2D tight-binding square lattice to be used 
as input for the various schemes. First and second moments of the irreducible partial densities 
of states are used to calculate these values up to cluster sizes of A^c = 8 {W = 4t is the half 
band- width). 



Nc 




Offset {iKjW) 


Width {ti/W) 


1 


(0,0) 


0.000 


0.500 


2 


(0,0) 


0.405 


0.293 




(vr, vr) 


-0.405 


0.293 


4 


(0,0) 


0.637 


0.218 




(^,0),(0,7r) 


0.000 


0.218 




(7r,7r) 


-0.637 


0.218 


8 


(0,0) 


0.811 


0.112 




(7r/2,7r/2),(37r/2,7r/2), 








(7r/2,37r/2),(37r/2,37r/2) 


0.000 


0.386 




(7r,7r) 


-0.811 


0.112 




(^,0),(0,7r) 


0.000 


0.181 



can display low dimensional properties in 2D. 

The essence of the proposed approach is shown for a cluster size of Nc = 4 in figure 
[2] for a quasi-2D system with in plane hopping, = 25tx (where is the inter-plane 
hopping). Panel (a) shows the numerically exact partial DOS for Nc = 4. Panel (b) shows 
the replacement of the partial DOS by Gaussians which match the first 2 moments. Panel 
(c) shows a modified finite size approximation, which replaces the partial DOS with a delta 
function matching the 1st moment only. The equivalent non-interacting "densities of states" 
associated with a finite size calculation are shown in panel (d). With a simple approximate 
form for the partial density of states, the cluster Green function can easily be calculated. I 
propose several approximate forms for the partial DOS for which an analytic form for the 
self-consistent condition can be found. These are summarized for convenience in table ITTl 
along with the analytic form of the Green function. A key feature of the approximation 
considered in this article is that it becomes better as cluster size increases since the second 
moment becomes vanishingly small, and the forms of the partial DOS converge. As such. 
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FIG. 2: Comparison between the partial non-interacting densities of states for (a) the DCA (cal- 



culated using the analytic tetrahedron method 



1) 



(b) Gaussian approximation to the DCA self- 



consistency (c) modified finite size calculations and (d) finite size calculations. A cluster size of 
Nc = 4 is shown since this is expected to be the most useful for analytic work. As the lattice size 
increases, the Gaussian form converges with the exact result. 

examinations of the four-site cluster {Nc = 4) represent the most challenging test of the 
approximate scheme. 

The Gaussian approximation to the partial DOS has the advantage that the full moment 
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TABLE II: Analytic forms for the partial DOS and the corresponding self-consistent equation, 



iu:n + IJ'—'^{Ki,i(jJn)-, s = sgn(Im[2:i]), erfc is the complimentary error function, ti = yc^. — {^kJ'^ 
and Cj is the average of the non-interacting dispersion over the course graining zone. 



Scheme 



partial DOS (A(e)) 



Green function {G{zi)) 



Gaussian 



exp(-(e-e-i)V2t?) 



tiV27r 



exp 



is(z-ei) - 



Semi-circular 



V4tf-(6-6-,)2 



2ntf 



{\ei -e\< 2ti 



2tf 



Square 



2tiV3 



In 



Z-€i+tiV3 



Lorentzian 



U 



^((e-e-i)Htf) 



1 



Zi-ei+itiSgn(lm[zi\) 



Modified FS 



S{€ - €i) 



1 

z—ei 



expansion is uniquely defined by the first two moments in the expansion, so all moments are 
well defined. However, the form of the Green function is relatively complicated, involving a 
complimentary error function. The semi-circular and square partial DOS are bounded, as is 
the case for all partial DOS in low dimensional systems, and lead to a Green function with 
a simple analytic form (i.e. the self-consistent equation can be inverted in terms of common 
functions). The first and second moments of the partial DOS are matched in all of these 
schemes. The square DOS approximation is not tested here. 

There are two additional schemes that have some potential for analytic calculations. The 
first is the Lorentzian approximation. For Nq — 1 this corresponds to the solution of an 
Anderson impurity model, and as such the self-consistent equations are decoupled. The 
second treats the partial DOS as a 5-function, where only the 1st moment is matched. This 
is similar to using a finite size (FS) technique, except that the shift in the peak location 
from the original position involves a large number of states, and such an approximation 
has aspects of both the thermodynamic limit of a large number of states and of finite size 
simulations. 



III. NUMERICAL TEST OF APPLICABILITY 



This section assesses the applicability of the approach described in the last section. The 
accuracy of the results computed using the approximate self-consistent step is considered 
as a function of the Hubbard U. Values of the Hubbard U approaching the band-width 
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FIG. 3: Comparison between the DCA and the approximate DCA schemes for thermodynamic 
quantities, using FLEX as a cluster solver and with a cluster size Nc = 4. Panel (a) shows the 
variation of the local moment {{Lq)) with coupling, and panel (b) shows the variation of the total 
energy, ^^totj with coupling. The temperature for these calculations is T = 0.04 and the density 
of electrons per site, n = 1. The Gaussian approximation to the DCA self-consistency leads to 
results which are almost indistinguishable from those calculated using the exact partial DOS (the 
open and filled circles lie on top of each other). The semi-circular approximation to the partial 
DOS also leads to good results. It is very important for the 1st and 2nd moments of the partial 
densities of states to match, or there are significant inconsistencies, such as those seen using the 
modified finite size approximation and Lorentzian DOS. For methods with matching moments, 
agreement is almost exact ior U ^W. 
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FIG. 4: Variation of filling with chemical potential showing the performance of the scheme away 
from half filling. T = 0.04H^ and U = 2W. Agreement is good for filling less than n = 1.5, where a 
kink develops. For larger fillings, there are small discrepancies of the order of a few percent. This 
is expected, since the Gaussian form of the partial DOS has unphysical long tails, where the exact 
partial DOS has a sharp cutoff. 



will be considered. In this paper, numerics are used to demonstrate the validity of the 
scheme. However, the intention is that researchers will be able to use approximate forms of 
the self-consistent condition to generate analytic solutions of model systems in conjunction 
with appropriate forms of the self energy. 

In order to test the accuracy of a scheme which uses approximate partial densities of 
states, the fluctuation-exchange (FLEX) approximation is used to compute the self energy 
of the cluster. Details of the FLEX approach and its use as a cluster solver may be found 



elsewhere 
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111 ]. FLEX has a transparent physical interpretation, since it approximates the 



self-energy by three subsets of diagrams, which represent spin-flips, density fluctuations and 
pair fluctuations. Computations were carried out at finite temperature using the Matsubara 
formalism. All calculations were initialized with zero self-energy, and then the Green function 
was calculated self-consistently until a convergence of 1 part in 10^ was reached. At this 
stage, observables were calculated. Ideally a full quantum Monte-Carlo scheme would be 
used to solve the cluster impurity problem, since FLEX can overestimate the self-energy 
and could potentially be insensitive to the different forms of the course graining schemes 
(for example the absence of van Hove singularities). However for moderate couplings FLEX 
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FIG. 5: Real and imaginary parts of the Green function on the Matsubara axis. Values are in 
agreement to within 1%. U = 0.5W, Nc = ^ and T = 0.04Ty and the system is half filled. 

is accurate, and it is expected that approximate schemes for the self-consistency are most 
likely to fail when electrons are not localized (since strongly localized electrons can be well 
described with a single site impurity). 

The total energy and local moment are plotted in figure [3] as coupling, f/, is varied. 
The cluster size is Nc = 4, temperature T = 0.04V1/^ and electron density n = 1. Results 
computed using a number of different approximations to the partial DOS are shown. The 
Gaussian approximation to the partial DOS leads to the most accurate results. In fact, the 
difference between the curves computed using the Gaussian and exact non-interacting DOS 
is not easily distinguished by eye. Agreement is also found when the semi-circular partial 
DOS is used to approximate the self-consistency, with an accuracy of a few percent. Using 
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FIG. 6: Comparison of real and imaginary parts of the Green function on the Matsubara axis. 
Parameters are as Fig. [5l 

the semi-circular DOS to approximate the self consistency is more favorable to analytic 
approximations, since the self-consistent equations also have a square root form and can 
easily be inverted. At large U, the curves corresponding to the Gaussian DOS, semi-circular 
DOS and the DCA scheme all converge. This is expected, since at large U at half-filling 
hopping is suppressed, the system becomes local and the precise form of the non-interacting 
DOS (which is related to the details of the intersite hopping and the specific lattice) becomes 
irrelevant. 

The matching of both moments is very important to achieve accurate results. When 
only the first moment is matched, for example using the modified finite size approximation, 
the errors on the moment and total energy are significant. There is a large overestimation 
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FIG. 7: Comparison of the density of states computed using a Fade approximant analytic contin- 
uation of order 800. While the detailed structure of the interacting DOS varies slightly, the main 
features are similar. This is in spite of the significant difference between the exact and approximate 
forms of the non-interacting density of states at the Fermi-surface. Farameters are as Fig. [5j 

of the local moment, and a huge underestimation of the total energy. Approximating the 
exact non-interacting density of states with a Lorentzian generates the worst approximation. 
Presumably this is because the second moment is ill defined, and the kinetic energy in the 
absence of interaction is grossly overestimated. However, some of the qualitative features 
remain and it might be possible to extract some useful physics even with such a crude 
approximation I2I. 

The previous results were calculated for a half filled lattice, corresponding to one of the 
interesting limits of the Hubbard model. To determine the performance of the approximation 
away from half-filling, figure H] shows the variation of the filling as the chemical potential, 
fi, is changed. The temperature is T = 0.04iy and the coupling U = 2W. Calculations 
are shown for the full DCA self consistency scheme, and also for the Gaussian and semi- 
circular approximations to the self-consistency. The change of n gives a rough measure of 
the variation in the renormalized DOS at the Fermi energy for a particular value of the 
chemical potential, and thus an indication of the performance of the scheme away from half 
filling. For n 7^ 1, it might be expected that the analytic schemes have larger errors, since 
the approximate forms of the non-interacting DOS do not have the correct form when the 
density of electrons is low (for example the Gaussian DOS has unphysical high energy tails). 
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In fact, the agreement is good for filling less than n = 1.5, where a kink develops. For larger 
fillings, there are only small discrepancies of the order of a fraction of a percent. 

Since it is possible that cancellations of errors could lead to good agreement of thermo- 
dynamic quantities, but may not give good results when correlation functions are computed. 
Fig. O shows the real and imaginary components of the Green function on the Matsubara 
axis. Only tiny differences can be seen between the forms of the Green functions computed 
using the Gaussian approximated scheme and the DCA. Here U = 0.51^, Nc = 4: and 
T = 0.04W^ and the system is half filled. To clarify the quality of the approximation, the 
difference between the values of the Green functions is shown in Fig. El For the param- 
eters used here, the difference is around 1% for the smallest Matsubara frequencies, and 
the two schemes converge for larger values of ujn- In the weak to intermediate U regime, 
any discrepancies between quantities computed using the approximate and exact forms of 
the self-consistency condition are expected to be most pronounced, and intermediate U is 
of most physical interest, so it is promising that the approximate scheme works well for 
U = 0.5W. 

Finally, to demonstrate that the form of the interacting density of states A{uj) = 
—lm[G{uj + ii])]/7i can be computed to reasonable accuracy using the alternative self con- 
sistency scheme, Fade approximants are used to analytically continue the Green function to 



the real axis 



12| as shown in Fig. [71 While the detailed structure of the interacting DOS 



varies slightly, the main features are similar. This is in spite of the significant difference 
between the exact and approximate forms of the non-interacting density of states at the 
Fermi-surface (for example, a factor of almost 2 difference can be seen at the Fermi surface 
in Fig. [2|). In particular, a similar amount of spectral weight is shifted across the Fermi 
surface in both cases. The density of states at the Fermi-surface associated with the course 
grained zone centered about the (vr, 0) point is of similar size for both approximations, and 
the band widening is similar. 



IV. EXTENSION TO TWO-PARTICLE PROPERTIES. 

In principle, the analytic schemes presented here could be used for calculating two par- 
ticle properties, such as the superconducting, charge density wave and spin density wave 
susceptibilities. To examine the potential for the computation of two particle properties. 
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I use the particle-hole susceptibility as an example. Following Ref. Il3|, the coarse-grained 
non-interacting susceptibility can be computed using the following expression, 

xl„,[{q,ius); {Ki,iuJn)iKj,iuJm) = Saa'Sij6nm^G„{Ki + k,iuJn)G„{Ki + k + q,iujn + ii^s) , 

k 

(5) 

where A; is a vector that lies within a sub-Brillouin zone. It is possible to simplify this 
expression for use with the approximate scheme discussed in this article in two special cases; 
when q = = (0, 0, ■ ■ ■ ) and when q = tt = (tt, vr, ■ ■ ■ ). First, if q = then the propagator 
becomes, 

E 

^ iun + - e(fc) - T^iiUn, Ki + k) iun + zz/s + yU - e(k) - T^Uun + ii^s, Ki + k) 

k 

Since in the dynamical cluster approximation S is constant within each sub-zone, this sim- 
phfies to. 



[iuJn + /i - e - Ki)){iuJn + ius + jJ. - e - T.{iuJn + iJ^s, Ki)) ' 

so that the form of the partial density of states consistent with the approximation to the 

self-consistent equation can be used. Equation [7] is appropriate for all forms of lattice. 

Simplifications can also be made in the more useful case where g = tt if there is particle- 
hole symmetry, for example if the lattice is square (or hypercubic) and only has near-neighbor 
hopping. Then e(fc) = —tJ2'i=i cos(/c), so e(fc + 7r) = — e(fc), and the susceptibility simplifies 

to 

xl^,[{7T,iiys);{Ki,iuj„){Kj,iuJm) = (8) 

nm J ^.^^ _^ ^ _ g _ Ki)){iuJn + iUs + jj, + e - + ius, Ki + tt)) ' 

The two particle cluster susceptibility Xc,aa'[{q,i^s)', {Ki,iLJn){Kj,iLJm)] is computed in 
a manner which is specific to the cluster solver that is used, and does not involve any 
coarse graining. Once the cluster properties have been computed, the coarse grained lattice 
susceptibility can be computed using, 

= x:' - ixT' + [xr\ (9) 



from which the spin and charge susceptibilities can be computed [l3|. Calculation of Xc is 
quite involved in the FLEX scheme and is not discussed here. 
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FIG. 8: Possible hopping scenarios for multiband models that can be treated using the approximate 
self-consistent scheme (a) two bands arise from sites of different energies in a lattice with a basis 
of two atoms (b) two bands arise from hopping between electrons in different orbitals. 

V. MULTIPLE BANDS 

I complete this article with a discussion of the applicability of the approximate scheme 
to multiple bands. I consider only two bands here, but similar arguments will be applicable 
when there are more than two bands. A very generic two-band Hamiltonian has the form, 

H2 = ^ tkcick + Ckdidk + 1k{c^kdk + (i^Cfc) (+ -ff2,int) , (10) 

k 

where and create electrons of distinct types (either on different atoms of a two-atom 
basis, or in different orbitals of a single atom), e^, is the electronic dispersion generated from 
hopping of electrons of type c, and Cfc corresponds to electrons of type d. The dispersions 
efc and Cfc typically correspond to hopping between orbitals of the same type on different 
sites as shown schematically as hops of amplitude t and t" in Fig. [Ht^b). 7fc is the dispersion 
of electrons that hop between different types of atom such as might be found in crystals 
with a two-atom basis (shown schematically in Fig. [H](a))- Dispersion of type 7^ can also 
originate from hopping between orbitals of different types on neighboring sites, as shown 
schematically as hops of amplitude t' in Fig. [8]^b). -?/2,int contains all the interaction terms 
between electrons of types c and d within the cell. 

Diagonalization of the quadratic terms using a standard Bogoliubov transformation 
yields, 

H2 = + Ef(3if3k (+i/2,int) , (11) 
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where, 



1 ± 



Cfc — Cfc 




2ii 



(12) 



V47^ + (efc - Cfc) 

4 = Ukai + VkPl d{ = -Vkai + UkPl ul = {l + Ak)/2, vl = (l-Afc)/2, u^Vk = y^l - Al/2 
and Ak = (efc — Cfc) / \/47| + (efc — Cfc)^- The non-interacting Matsubara Green functions for 
r > form a 2 x 2 matrix with the elements, G^cc{k,T) = — (cfc(r)c|,(0)) and GfJ{k,T) = 
— (cfc(T)c/^(0)) etc. which can be computed by substituting the Bogohubov transformation 
(since a and /3 excitations have a Fermi-Dirac distribution). After Fourier transformation. 



and 



iujr, — E, 



UkVk 



(a) 



+ 



iUJ„ - E; 



(/3) 



E 



(a) 



E 



(/3) 



(z^„-4"))(z^.-Ef )■ 
When interaction is switched on, the coarse grained Green function is. 



(13) 



(14) 



(15) 



where the self-energy, S, is also a 2 x 2 matrix. As in the one-band case, the self-energy is 
constant within each sub-zone. 

In contrast to the one-band case, the sum on k cannot always be replaced by a density 
of states. This is because e^, Cfc and jk do not normally have the same form. Under certain 
conditions it is possible to make the replacement and therefore use the analytic forms of 
the self-consistency considered in this article. In the first case, bands with dispersion 7^, are 
formed via nearest-neighbor hops in a system with a basis of two atoms, and and Cfc are 
flat, representing the energies of the atomic sites. Then the sum over momenta is replaced 
by an integral over the density of states corresponding to 7fc, which can then be replaced 
using one of the simplified analytic forms considered in section [III I note that two-band 
systems of this type have been studied using DMFT Q- 

The sum over momenta can also be replaced with an integral over a density of states 
when Cfc oc Cfc, Cfc •jk or Cfc oc ■jk (up to an offset) so long as any other dispersions are 
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flat. This is possible for electrons moving between atoms with several orbitals by nearest 
neighbor hopping only (the hopping integrals can be different) as shown schematically in Fig. 
MJo)- Then the density of states corresponding to only one of the dispersions can be used, 
since any dispersion can be represented in terms of the other dispersion. In general systems 
with next-nearest neighbor hopping are excluded from this argument unless there are good 
grounds for expecting the ratio of nearest to next-nearest neighbor hopping integrals to be 
identical for every band. There may be other more complicated non-linear maps between 
dispersions (e.g. = -F(7fc) where F is an arbitrary function) so that the same trick can be 
applied, however the mapping must be k independent so such mappings will also be special 
cases. 



VI. SUMMARY 



This paper introduced approximate schemes for the self- consistent step of the dynamical 
cluster approximation. In order to simplify the self-consistent step, the partial density of 
states used to calculate the Green function was replaced with a simple analytic form. The 
limitations and errors associated with the approach were tested using the FLEX approxi- 
mation for the self energy. The approximate schemes were demonstrated to work well for 
intermediate coupling over a range of fillings. 

The best approximate result comes from replacing the exact non-interacting DOS with 
a Gaussian. Approximating using a semi-circular partial DOS also leads to agreement of 
thermodynamic properties to within a few percent of the exact result. The semi-circular ap- 
proximation leads to self-consistent equations with a much simpler form, which is promising 
for analytic work. Finally, forms such a Lorentzian and a modified finite size scheme, where 
only the 1st moment is matched, lead to quantitatively different results, although some 
of the qualitative physics remains. These schemes lead to even simpler forms for the self- 
consistency. Examination of the Green function and carrying out an analytic continuation 
demonstrate that the scheme can also compute dynamic properties to reasonable accuracy. I 
add the caveat that some properties of the 2D Hubbard model may show features consistent 
with the influence of van Hove singularities and such singularities are not present in 
the simple approximations with small cluster sizes (although they will emerge in very large 
clusters). In 3D, there are no divergences in the non-interacting density of states and spatial 
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fluctuations are less relevant, so the schemes may be most effective for such problems. 

Given that the analytic forms for the self-consistency introduced here work effectively, 
the formalism could provide a good starting point for the calculation of analytic phase 
diagrams and other properties in low dimensional systems where non-local fluctuations are 
essential (for example unconventional superconductivity such as d, p and extended-s wave 
which cannot be examined with the DMFT) . It is hoped that this article will stimulate such 
studies. 
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